function [ contribERdC ] = ERdC_contrib_M2( donnees,elementCourant )

E = donnees.poutre.E ;
I = donnees.poutre.I ;
L = donnees.poutre.L ;
F = donnees.F ;

a = donnees.forceLineiqueA ;
b = donnees.forceLineiqueB ;
c = donnees.forceLineiqueC ;

xi = elementCourant(1) ;
xj = elementCourant(2) ;

%intM2 = (1/144)*(-((6/5)*xi^5*(L*(a^2*L^3 - 4*a*(b*L^2 + F) - 8*b^2*L) + 2*c*L*(a*L - 8*b) + 16*b*F + 6*c^2)) + (6/5)*xj^5*(L*(a^2*L^3 - 4*a*(b*L^2 + F) - 8*b^2*L) + 2*c*L*(a*L - 8*b) + 16*b*F + 6*c^2) - (a^2*xi^9)/9 + (a^2*xj^9)/9 - (4/7)*xi^7*(3*a*c + 4*b^2) + (4/7)*xj^7*(3*a*c + 4*b^2) - (4/3)*xi^3*(L^2*(3*c*L*(11*a*L + 16*b) + L^2*(2*a*L + 3*b)^2 + 54*c^2) - 12*F*L*(L*(2*a*L + 3*b) + 9*c) + 36*F^2) + (4/3)*xj^3*(L^2*(3*c*L*(11*a*L + 16*b) + L^2*(2*a*L + 3*b)^2 + 54*c^2) - 12*F*L*(L*(2*a*L + 3*b) + 9*c) + 36*F^2) + 2*xi^4*(3*c*(L^2*(2*a*L + b) - 6*F) + b*L*(12*F - L^2*(3*a*L + 4*b)) + 18*c^2*L) - 2*xj^4*(3*c*(L^2*(2*a*L + b) - 6*F) + b*L*(12*F - L^2*(3*a*L + 4*b)) + 18*c^2*L) - L^2*xi*(L*(L*(3*a*L + 4*b) + 6*c) - 12*F)^2 + L^2*xj*(L*(L*(3*a*L + 4*b) + 6*c) - 12*F)^2 - (2/3)*xi^6*(3*b*(4*c - a*L^2) - 2*a*(a*L^3 + 3*c*L - 3*F)) + (2/3)*xj^6*(3*b*(4*c - a*L^2) - 2*a*(a*L^3 + 3*c*L - 3*F)) + 2*L*xi^2*(L*(L*(2*a*L + 3*b) + 6*c) - 6*F)*(L*(L*(3*a*L + 4*b) + 6*c) - 12*F) - 2*L*xj^2*(L*(L*(2*a*L + 3*b) + 6*c) - 6*F)*(L*(L*(3*a*L + 4*b) + 6*c) - 12*F) - a*b*xi^8 + a*b*xj^8) ;

%intM2 = (1/144)*((-L^2)*(12*F + L*(6*c + L*(4*b + 3*a*L)))^2*xi + 2*L*(6*F + L*(6*c + L*(3*b + 2*a*L)))*(12*F + L*(6*c + L*(4*b + 3*a*L)))*xi^2 - (4/3)*(36*F^2 + 12*F*L*(9*c + L*(3*b + 2*a*L)) + L^2*(54*c^2 + L^2*(3*b + 2*a*L)^2 + 3*c*L*(16*b + 11*a*L)))*xi^3 + 2*(18*c^2*L + 3*c*(6*F + L^2*(b + 2*a*L)) - b*L*(12*F + L^2*(4*b + 3*a*L)))*xi^4 - (6/5)*(6*c^2 - 8*b^2*L^2 + 2*c*L*(-8*b + a*L) - 4*b*(4*F + a*L^3) + a*L*(4*F + a*L^3))*xi^5 + (2/3)*(3*b*(-4*c + a*L^2) + 2*a*(3*F + 3*c*L + a*L^3))*xi^6 - (4/7)*(4*b^2 + 3*a*c)*xi^7 - a*b*xi^8 - (a^2*xi^9)/9 + L^2*(12*F + L*(6*c + L*(4*b + 3*a*L)))^2*xj - 2*L*(6*F + L*(6*c + L*(3*b + 2*a*L)))*(12*F + L*(6*c + L*(4*b + 3*a*L)))*xj^2 + (4/3)*(36*F^2 + 12*F*L*(9*c + L*(3*b + 2*a*L)) + L^2*(54*c^2 + L^2*(3*b + 2*a*L)^2 + 3*c*L*(16*b + 11*a*L)))*xj^3 - 2*(18*c^2*L + 3*c*(6*F + L^2*(b + 2*a*L)) - b*L*(12*F + L^2*(4*b + 3*a*L)))*xj^4 + (6/5)*(6*c^2 - 8*b^2*L^2 + 2*c*L*(-8*b + a*L) - 4*b*(4*F + a*L^3) + a*L*(4*F + a*L^3))*xj^5 - (2/3)*(3*b*(-4*c + a*L^2) + 2*a*(3*F + 3*c*L + a*L^3))*xj^6 + (4/7)*(4*b^2 + 3*a*c)*xj^7 + a*b*xj^8 + (a^2*xj^9)/9) ;

intM2 = (1/144)*((-L^2)*(12*F + L*(6*c + L*(4*b + 3*a*L)))^2*xi + 2*L*(6*F + L*(6*c + L*(3*b + 2*a*L)))*(12*F + L*(6*c + L*(4*b + 3*a*L)))*xi^2 - ...
   (4/3)*(36*F^2 + 12*F*L*(9*c + L*(3*b + 2*a*L)) + L^2*(54*c^2 + L^2*(3*b + 2*a*L)^2 + 3*c*L*(16*b + 11*a*L)))*xi^3 - (-36*c^2*L - 12*c*(3*F + L^2*(b + a*L)) + b*L*(12*F + L^2*(4*b + 3*a*L)))*xi^4 - ...
   (2/5)*(18*c^2 - 12*b^2*L^2 + 6*c*L*(-4*b + a*L) + 3*a*L*(4*F + a*L^3) - 4*b*(6*F + a*L^3))*xi^5 - (2/3)*(b*(6*c - 3*a*L^2) - 2*a*(3*F + 3*c*L + a*L^3))*xi^6 - (4/7)*(b^2 + 3*a*c)*xi^7 - (1/2)*a*b*xi^8 - (a^2*xi^9)/9 + ...
   L^2*(12*F + L*(6*c + L*(4*b + 3*a*L)))^2*xj - 2*L*(6*F + L*(6*c + L*(3*b + 2*a*L)))*(12*F + L*(6*c + L*(4*b + 3*a*L)))*xj^2 + ...
   (4/3)*(36*F^2 + 12*F*L*(9*c + L*(3*b + 2*a*L)) + L^2*(54*c^2 + L^2*(3*b + 2*a*L)^2 + 3*c*L*(16*b + 11*a*L)))*xj^3 + (-36*c^2*L - 12*c*(3*F + L^2*(b + a*L)) + b*L*(12*F + L^2*(4*b + 3*a*L)))*xj^4 + ...
   (2/5)*(18*c^2 - 12*b^2*L^2 + 6*c*L*(-4*b + a*L) + 3*a*L*(4*F + a*L^3) - 4*b*(6*F + a*L^3))*xj^5 + (2/3)*(b*(6*c - 3*a*L^2) - 2*a*(3*F + 3*c*L + a*L^3))*xj^6 + (4/7)*(b^2 + 3*a*c)*xj^7 + (1/2)*a*b*xj^8 + (a^2*xj^9)/9) ;

%DeltaM2 = intM2 - (-(1/3))*F^2*(xi - xj)*(3*L^2 + xi^2 + xi*xj + xj^2 - 3*L*(xi + xj))

contribERdC = 1/(E*I) * intM2 ;

end

